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ABSTRACT.  This  report  describes  the  theory  and 
design  of  an  optical  ozonesonde  system  developed 
for  use  on  Areas  rockets.  The  payload  is  an 
ultraviolet  filter  photometer  useful  in  the  alti¬ 
tude  range  20  to  60  km.  A  real-time  digital  data 
logging  system  has  been  constructed  for  recording 
the  data  that  are  telemetered  on  the  1680  mHz 
meteorological  frequency.  Data-processing  proce¬ 
dures  are  described  and  results  of  a  test  flight  J' 
are  presented. 
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FOREWORD 

This  report  describes  the  design,  development,  and  testing  of  a 
system  for  measurement  of  the  distribution  of  atmospheric  ozone  using 
meteorological  rockets. 

The  objective  of  this  program  was  to  produce  a  light,  compact, 
rugged,  and  inexpensive  ozonesonde  (ROCOZ)  for  use  with  Areas  rockets 
to  collect  data  in  the  region  from  20  to  60  km.  Such  an  instrument 
was  required  to  determine  if  ozone  is  in  photocheaiical  equilibrium 
with  oxygen  and  ultraviolet  sunlight  as  the  present  theory  predicts. 

This  work  is  part  of  a  program  directed  to  the  measurement  of 
physical  characteristics  and  composition  of  the  stratosphere  to  an 
altitude  of  60  km.  The  work  complements  the  intensive  balloon 
investigations  of  lower  stratospheric  ozone  and  bears  a  direct 
relationship  to  the  Meteorological  Rocket  Network  measurements  of  air 
temperature  and  vinds. 

Operational  flights  have  been  conducted  at  geographical  locations 
of  interest  and  results  will  be  reported  ir  future  publications. 

The  project  was  initially  supported  by  the  Office  of  Naval  Rasaarch 
under  Project  Number  NR  082-183  and  then  by  the  Pacific  Missile  Range 
Project  Orders  PO-3-0021/072831005  and  PO-4-0005/07 3146202.  Rockets 
and  parachutes  were  furnished  by  the  Geophysics  Division  of  PMR. 
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INTRODUCTION 
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The  rocket  ozonesonde  (ROCOZ),  Fig.  1,  developed  at  this  Center  is 
an  optical  instrument  using  the  sun  as  a  light  source  and  four  selected 
pass  bands  in  the  ultraviolet  region  from  2700  to  3300  A  to  measure 
attenuation  of  solar  energy  by  ozone.  An  Areas  rocket  carries  the  ROCOZ 
to  an  altitude  of  55  to  65  km  and  a  gas  generator  ejects  the  payload  for 
deployment  on  a  parachute.  The  relatively  slow  rate  of  descent  allows  a 
longer  time  for  measurement  than  vould  be  obtained  in  free  fall. 

ROCOZ  is  designed  as  a  sequential  sampling  filter  photometer  that 
makes  use  of  a  standard  meteorological  frequency  and  ground  station  for 
telemetry.  An  automatic  digital  recording  system  was  constructed  to 
receive  the  data  on  the  ground. 

The  basic  system  operated  successfully  in  test  flights.  However, 
the  tests  did  show  that  the  payload-parachute  system  must  be  stabilized 
to  acquire  higher  quality  analog  data  from  this  payload  configuration 
over  the  entire  altitude  range. 


THEORETICAL  CONSIDERATIONS 


OPTICAL  TECHNIQUE  FOR  ATMOSPHERIC  OZONE  MEASUREMENTS 

Ozone  is  characterized  by  an  intense  optical  absorption  band  in  the 
near  ultraviolet.  As  measured  by  Inn  and  Tanaka  (Ref.  1),  the  exponen¬ 
tial  absorption  coefficient  reaches  a  peak  value  of  311  cm-*  (base  e) 
at  2553  A,  but  decreases  to  less  than  1  cm"1  for  wavelengths  greater 
than  3180  A.  With  a  typical  total  ozone  amount  of  0.3  cm  n.t.p., 
essentially  no  radiation  reaches  the  earth's  .surface  at  wavelengths 
less  than  2950  A.  The  total  ozone  amount  in  a  column  through  the 
atmosphere  is  measured  routinely  at  more  than  10C  stations  around  the 
world  with  the  Dobson  spectrophotometer  and  similar  instruments  by 
optical  absorption  at  wavelengths  longer  than  3050  A.  These  instruments 
compare  the  irradiance  from  extraterrestrial  sources  (primarily  the  sun) 
at  wavelength  pairs  separated  by  approximately  200  A. 
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A  similar  approach  has  been  used  to  determine  the  height  distribution 
of  ozone  with  balloon-borne  spectrographs  and  photometers.  Paetzold  in 
Germany,  Vassy  in  France,  and  Kobayashi,  et  al.,  in  Japan  have  made  numer¬ 
ous  soundings  with  filter  photometers  (Ref.  2-4).  The  Paetzold  instrument 
is  equipped  to  telemeter  data  on  solar  irradlance  in  two  spectral  bands 
between  3000  and  4000  X  to  ground  stations.  One  band  is  located  in  the 
ozone  absorption  region;  the  other  band  where  absorption  is  zero.  The 
ozone  concentration  as  a  function  of  height  is  calculated  from  the  change 
irr-’diance  in  the  absorbed  band  with  height.  To  avoid  having  to  point 
the  optical  axis  these  instruments  are  constructed  with  a  wide  field  of 
view.  This  results  in  a  moderation  of  the  irradlance  gradient  at  alti¬ 
tudes  where  the  skylight  is  a  significant  component  of  the  total  irra- 
d lance.  This  is  generally  the  case  below  15-20  km,  depending  on  the 
separation  between  the  bands.  The  problem  is,  of  course,  most  difficult 
at  the  ground  where  the  Dobson  spectrophotometer  is  used.  Here  both 
Rayleigh  and  Mie  scattering  provide  large  contributions  to  the  skylights 
Only  by  restricting  the  field  of  view  and  comparing  overlapping  wave¬ 
length  pairs  is  the  problem  minimized. 

From  20  km  to  normal  maximum  balloon  altitudes  the  optical  ozoneeondes 
are  capable  of  accurate  measurement  of  the  broad  details  of  the  ozone 
distribution.  The  ability  to  determine  detailed  structure  is  related  to 
the  accuracy  of  the  height  gradient  of  irradlance.  With  typical  photom¬ 
eter  and  telemetry  noise  levels,  data  over  at  least  1  km  are  required  to 
establish  the  gradient,  and  the  resolution  is  typically  1  to  2  km,  which 
for  most  research  purposes  is  adequate. 


SELECTION  OF  WAVELENGTHS  FOR  HIGH  ALTITUDE  MEASUREMENTS 

In  the  balloon  altitude  regime,  x(h),  the  total  ozone  remaining  above 
height,  h,  ranges  from  a  normal  r— <imura  of  0.5  cm  at  the  ground  to  0.1  cm 
at  30  kn.  The  solar  irradlance,  j.(A,h),  at  yavelength,  X,  and  height,  h, 
depends  on  the  absorption  coefficient  i(A)  and  x(h),  acccrd_ng  to  Beer's 
law: 

J(A,h)  ■  I(a,«)  exp[-a(A)  x(h)  sec  z] 

where  z  is  the  solar  zenith  distance  and  I (A,*)  is  the  extraterrestrial 
solar  irradlance.  For  low  values  of  z,  the  relative  irradlance  at  3100  X 
varies  from  0.25  at  the  ground  to  0.75  at  30  km.  Above  the  primary  ozone 
peak  at  25  km  the  concentration  is  predicted  by  photochemical  theory  to 
decrease  approximately  exponentially  with  height.  Early  measurements  on 
rockets  using  spectrographic  payloads  have  indicated  that  this  is  correct. 
Whereas  a  single  filter  at  3100  X  (and  a  reference  at  a  longer  wavelength) 
has  been  adequate  for  balloon  ozonesondes,  this  is  no  longer  the  case  at 
higher  altitudes  where  x(h)  is  expected  to  decrease  by  three  orders  of 
magnitude  between  30  and  60  km. 
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If  it  is  assumed  that  the  ozone  concentration,  c(h),  above  the 
primary  peak  can  be  represented  by  the  function 

e(h)  -  tQ  exp ( - (h  -  hQ)/H] 

where  H  is  the  ozone  scale  height,  then  the  total  ozone  (above  height  h) 
is  given  by 

x(h)  -  cQ  H  exp[-(h  -  hQ)/H]  (2) 

The  ozone  concentration  is  determined  from  the  change  in  irradiance  with 
height.  From  Eq.  1  and  2,  the  following  can  be  written  for  any 
wavelength 

3 1 

7T*  a  t  I  sec  z  exp[-(h  -  h  )/H] 
on  o  o 

This  function  has  a  maximum  at  h  *  h  Hi  In  U  a  c  sec  z. 

mo  o 

If  h^  is  selected  as  the  reference  height,  it  is  found  that 

-p  *  exp[-u  -  exp(-u) J 
o 

where  u  *  (h  -  h^j/H.  This  function,  shown  in  Fig.  2,  was  originally 
derived  by  Chapman,  Ref.  5,  to  show  that  optical  absorption  In  a  plane 
exponential  atmosphere  takes  place  in  relatively  thin  strata  when  the 
absorption  coefficient  is  constant.  In  the  present  context,  the  maxi¬ 
mum  useful  height  range  for  a  narrow-band  optical  filter  can  be  deter¬ 
mined  from  this  curve.  For  example,  assuming  a  scale  height  of  6  km 
and  r.  minimum  detectable  change  in  irradiance  of  2Z/km,  3l/3h  H/I  *  0.12. 
The  absorption  rate  function  exceeds  this  value  between  u  *  -1.24  and 
+1.98.  This  corresponds  to  a  total  height  range  of  19.3  km.  These 
assumptions  are  not  unrealistic  and  it  can  be  expected  that  *_.,e  ir radi¬ 
ance  at  two  wavelengths  corresponding  to  different  valves  of  h  will  be 
required  for  ozone  determinations  over  the  altitude  range  from  30  to 
60  km  where  the  concentration  decays  exponentially  with  height.  The 
optimum  values  of  h  are  in  the  vicinity  of  33  and  47  km.  The  wave- 

III 

lengths  required  to  obtain  these  optimum  heights  can  be  calculated  given 

the  a (a)  curve  and  an  estimate  of  x(h)  sec  z  from  the  condition: 

a  (A)  x(h  )  sec  z  ■  1. 

01 
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Assuming  a  nominal  z  of  45  degrees 

a(Xj)  ■  l/x(33  km)  sec  45*  ■  1/.05  x  1.414  ■  14.1  cm-1 

a(X2)  -  l/x(47  km)  sec  45*  -  1/.003  x  1.414  -  238  cm-1 

Following  recommendations  of  Inn  and  Tanaka  (Ref.  1),  their  measurements 
of  a  (A)  have  been  adopted  for  wavelengths  less  than  2750  A  and  the 
measurements  of  Vigroux  (Ref.  6)  have  been  adopted  for  longer  wavelengths 
because  of  the  higher  resolution.  The  absorption  coefficients  derived 
above  correspond  to  center  wavelengths  of  2970  A  at  33  km  and  2675  A 
at  47  km.  While  these  are  optimum  wavelengths  some  adjustment  is  made 
from  consideration  of  the  solar  spectrum  ana  the  actual  height  range 
attainable  with  the  Areas  system. 

For  the  region  below  30  km  a  band  near  3100  A  is  suitable  as  used  in 
the  Paetzold  balloon  sonde.  For  a  reference  filter  where  a (A)  -  0  a 
center  wavelength  near  3300  A  has  been  adopted. 


CALCULATION  OF  OZONE  CONCENTRATIONS  FROM  THE  OBSERVED  IRRADIANCE  DATA 


In  the  calculation  of  the  ozone  distribution  from  irradiance  data. 
Beer's  law  (Eq.  1)  is  used.  However,  the  optical  filters  have  bandwidths 
wide  enough  that  the  wavelength  dependence  of  various  input  parameters 
must  be  recognized.  It  is  possible  to  treat  the  ozone  absorption  coeffi¬ 
cient  and  Rayleigh  extinction  coefficient  as  functions  of  the  total  ozone 
amount,  x,  and  the  relative  air  mass*  m,  because  these  are  the  parameters 
that  result  in  an  alteration  of  the  spectrum  as  a  function  of  height.  An 
"effective  absorption  coefficient"  is  defined  by: 


<Xj(x,m)  -  -  -  In 


A2 

Jx  I (A,®)  Fj(A)  P(A)  t (A)  exp [-a (A)  x-B(A)m]dA 


T 


(3) 


I(A,«)  Fj  (A)  P(A)  t (A)  exp(~8(A)m]dA 


where 

til 

F^(a)  *  transmission  function  of  the  j  filter 

t(X)  *  transmission  function  of  the  integrating  sphere 
P(A)  *  photodetector  response 
0 (A)  ■  Rayleigh  extinction  coefficient 


and  the  wavelength  limits,  Aj  and  cover  the  wavelength  region  of 
interest.  Because  of  variations  between  filters  constructed  for  the 
various  payloads,  it  has  been  necessary  to  calculate  for  each  new 
filter.  Typical  results  are  shown  in  Fig.  3. 
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FIG.  3.  Absorption  Coefficients  of  Filters. 


The  basic  absorption  equation  can  then  be  written  as  the  ratio  of 

the  output  signals,  S.(h),  for  the  jth  filter  and  the  reference  filter 

(j  -  o).  J 


Sj  (h) 

so(h) 


s  (-) 
s  (®) 

O 


expy-[cXj  (x,m)  -  oo(x,m)]  x(h)  sec  z 
exp{-[ej(x,m)  -  Bo(x,m)]  ra(h)| 


(4) 


where  Bj(x,m)  is  defined  similarly  to  Oj(x,m).  The  extraterrestrial 
ratio  (h  ■  ®)  can  be  measured  on  a  rocket  flight  in  that  the  value  of 
ajx  =  0  at  the  top  of  the  flight.  At  low  altitudes  (h  <  20  lan)  sky¬ 
light  begins  to  contribute  to  the  total  irradiance  as  observed  by 
Paetzold  (Ref.  2).  In  the  present  analysis  where  the  interest  is 
primarily  in  the  ozone  concentrations  in  and  above  the  primary  peak 
this  contribution  can  be  ignored  for  the  present  time.  As  data  are 
accumulated  with  the  filters  selected  for  this  instrument,  attempts 
will  be  made  to  derive  corrections  for  application  to  the  lower 
altitude  data. 
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The  ozone  concentrations,  Ax(h)/Ah,  are  calculated  using  a 
simplified  form  of  a  finite  difference  equation  as  derived  by  Paetzold 
(Ref.  2): 


Ax(h) 

Ah 


-  a  )  sec  z  “ 
o 


A  In  S,(h)/S  (h) 

_ J  ° 

Ah 


+ 


Ad(h,z)  8aj 
Ah  3d  (h) 


d{h,z) 


A  sec  z 
Ah 


+ 


Am(h,z)  3aj 
Ah  3m(h,z) 


(5) 


where  d  *  x(h)  sec  z.  This  equation  yields  a  solution  by  iteration 
given  an  assumed  starting  distribution,  dj(h,z).  The  equation  is 
greatly  simplified  at  high  altitudes  where  m  =  0.  Also  in  parachute 
descent  the  altitudes  from  60  to  30  km  are  covered  in  approximately 
10  min  so  that  sec  z  is  essentially  constant  for  low  z. 


Equation  5  is  strictly  true  only  if  the  ozone  concentration  is 
approximate  /  constant  with  height.  Above  the  primary  ozone  peak  the 
concentration  decays  roughly  exponentially  with  height  and  a  correction 
can  be  applied  to  obtain  the  true  concentration,  e(h).  In  this  case, 


e(h  ) 
o 


Ax/ Ah 


H[exp(-  f)  -  l' 


(6) 


where  h  is  the  base  of  the  altitude  interval.  Ah,  and  H  is  the  ozone 
scale  height  (assumed  constant  over  Ah) . 


R0C0Z  SYSTEM  DESIGN 


FLIGHT  CHARACTERISTICS 

Atmospheric  soundings  using  the  Areas  rocket  generally  use  a 
dropsonde  technique.  The  rocket  carries  a  payload  to  an  altitude  of 
55  to  65  km  where  a  gas  generator  ejects  the  payload  for  deployment  on 
a  parachute.  As  mentioned  earlier,  the  relatively  slow  rate  of  descent 
allows  a  longer  time  for  measurement  than  would  be  obtained  in  free 
fall.  (A  typical  height  vs.  time  profile  is  shown  on  page  26.)  The 
rocket  is  launched  from  a  closed-breech  tube  that  creates  an  initial 
speed  of  150  ft/sec  to  decrease  wind  effects  on  the  motor-case  impact 
location.  A  200-msec  launch  shock  peaking  near  60  &  is  found.  The 
motor  burns  for  28  sec  giving  the  system  an  acceleration  of  A  to  5  £. 

At  burnout  a  100-sec  pyrotechnic  delay  train  is  initiated  that,  in  turn, 
ignites  the  ejection  mechanism  near  apogee.  The  ejection  system  produces 
a  brief  70-£  shock. 
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A  15  ft  metallized  silk  parachute  has  generally  been  used  with  this 
system.  The  descent  rate  is  proportional  to  air  density,  varying  from 
135  m/sec  at  55  km  to  20  m/sec  at  30  km.  During  this  portion  of  descent 
a  high  amplitude  pendulation  is  observed.  This  is  due  to  an  Initial 
1,200-rpm  spin  required  for  rocket  stability  and  the  parachute  character¬ 
istics.  The  pendulation  is  gradually  damped  out  by  aerodynamic  drag  and 
internal  friction  in  the  parachute.  Below  30  km  pendulation  due  to  wind 
shears  is  observed. 

Scientific  data  from  the  payload  are  generally  transmitted  on  a 
1680-mHz  carrier  to  a  GMD  meteorological  receiver.  Information  on  alti¬ 
tude  and  winds  is  obtained  by  a  tracking  radar. 


DESIGN  CRITERIA 

The  basic  task  is  to  design  a  multichannel  ultraviolet  filter 
photometer,  capable  of  stable  operation  in  the  environmental  conditions 
set  by  the  rocket  flight  followed  by  slow  parachute  descent  through  the 
atmosphere.  In  addition  to  mechanical  conditions  stated  above,  the 
payload  must  incorporate  thermal  insulation  as  it  is  exposed  to  atmos¬ 
pheric  temperatures  ranging  from  +30*  C  to  -80*C.  The  package  design  is 
thus  somewhat  restricted. 

A  second  criterion  has  been  the  cost  of  expendables.  The  general 
rule  has  been  that  the  payload  cost  should  not  exceed  the  vehicle  cost 
by  a  large  margin.  Because  system-use  costs  also  involve  data  readout 
and  reduction,  a  great  advantage  is  gained  by  automatic  digitizing  of 
the  data.  For  this  reason  the  system  has  been  designed  as  a  point 
sampler,  incorporating  a  signal  that  allows  direct  real-time  digitization 
of  the  signal.  The  optical  sampling  rate  is  determined  by  the  desired 
point  vs.  height  density  and  the  maximum  descent  rate.  Given  a  minimum 
point  density  of  5  per  kilometer,  one  sample  per  second  of  each  optical 
channel  is  required  from  the  flight  profile. 

A  third  criterion  is  field  serviceability.  The  payload  is  designed 
for  convenient  disassembly  and  reassembly  for  replacement  of  modular 
components  while  maintaining  optical  alignment. 


FLIGHT  INSTRUMENT  DESIGN 

The  ozonesonde  is  designed  as  a  sequential  sampling  filter 
photometer.  It  makes  use  of  a  standard  meteorological  frequency  and 
ground  station  for  telemetry  but  requires  specialized  equipment  for  data 
recording  on  the  ground.  To  provide  flexibility  in  channel  capacity  and 
bandwidth,  subcarrier  oscillators  and  CW  modulation  have  been  utilized. 
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A  block  dlagran  of  Che  sounding  payload  Is  shown  in  Fig.  A.  In  order 
to  obtain  the  required  data  point  density  consistent  with  data  readout 
capability,  an  automatic  digital  recording  system  has  been  constructed 
with  real-time  capability. 


FIG.  A.  Block  Diagram  for  ROCOZ  Gzonesonde. 


PHOTOMETER  CONSTRUCTION 

Integrating  Sphere.  To  avoid  the  complexity  of  an  active  pointing 
system  the  photometer  has  been  designed  to  accept  light  from  a  wide 
angle  field  of  view.  The  optical  system  is  adapted  from  that  used  in 
the  Paetzold  ozonesonde.  As  in  that  instrument  a  transmission-type  inte¬ 
grating  sphere  accepts  light  from  all  unobscured  angles.  The  device 
consists  of  a  hollow  Suprasll  sphere,  open  on  the  base,  and  internally 
coated  with  an  optically  thick  layer  of  magnesium  oxide.  Light  incident 
on  the  surface  of  the  spherical  portion  is  scattered  by  the  magnesium 
oxide  coating.  The  portion  that  is  scattered  into  this  sphere  is 
internally  scattered  and  finds  its  way  to  the  photometer  entrance 
aperture.  The  light  available  for  the  photometer  is  dependent  geometri¬ 
cally  on  the  illuminated  surface  area.  For  a  sphere  this  is  independent 
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of  the  angle  of  incidence  so  the  performance  is  controlled  by  the 
quality  of  the  magnesium  oxide  coating.  A  properly  constructed  sphere 
has  a  homogeneous  coating  with  a  direct  ray  transmission  of  0.01  to 
0.05X  and  a  net  transmission  of  1Z  to  2Z. 

With  a  semi-isotropic  receiver,  at  some  altitude  the  scattered  light 
becomes  a  measurable  fraction  of  the  total  light  and  a  correction  becomes 
necessary  depending  on  the  wavelength.  Because  most  of  this  light  comes 
from  angles  near  and  below  the  horizon,  a  shield  limiting  the  field  of 
view  to  elevation  angles  greater  than  either  20  or  30  degrees  has  been 
installed.  This  has  a  second  function  of  protecting  the  sphere  during 
the  payload  ejection  sequence. 

Optical  Filters.  The  optical  filters  are  physically  divided  into 
two  sections:  (1)  a  broad-band  filter  that  serves  to  Isolate  the  ultra¬ 
violet  region  of  interest  and  effectively  block  the  visible  region,  and 
(2)  a  set  of  four  narrow-band  filters  that  define  the  particular  bands 
of  interest.  The  broad-band  filter  consists  of  12-mn  NiSO^'b^O  sand¬ 
wiched  between  3-mm  plates  of  Corning  9863  glass.  An  absorbance  curve 
of  a  representative  filter  is  shown  in  Fig.  5. 

The  narrow-band  filters  are  composite  interference  and  chemical 
filter  assemblies  constructed  to^have  half-widths  of  40-50  A  centered 
near  2700,  2950,  3050,  and  3300  A.  Transmission  curves  for  a  repre¬ 
sentative  set  are  shown  in  Fig.  6.  At  the  present  time  these  filter 
sets  are  constructed  in-house  using  commercial  components  where  avail¬ 
able.  Primary  criteria  are  center  wavelength,  half-width,  and  equiva¬ 
lent  band  pass  or  peak  transmission.  The  latter  parameter  is  determined 
such  as  to  provide  equal  signals  from  all  filters  outside  the  ozone 
region  (h  >  60  km). 

The  narrow-band  filters  are  installed  in  a  five-position  filter 
wheel  in  which  one  position  is  filled  with  an  opaque  insert.  This 
wheel  is  rotated  continuously  at  60  rpm  by  a  Brailsford  motor  mounted 
inside  a  magnetic  shield  to  minimize  radiation.  The  optical  components 
are  shown  in  Fig.  7. 

Photodetector  Selection.  The  types  of  photodetectors  available  for 
use  in  the  near  ultraviolet  are  unfortunately  veiy  limited.  Photomulti¬ 
pliers  and  vacuum  photodiodes  are  bulky  and  expensive  but  are  quite 
desirable  in  terms  of  dynamic  range  and  linearity.  Of  the  solid-state 
photodetectors  only  special  selenium  photocells  and  some  varieties  of 
silicon  cells  have  residual  sensitivity  down  to  2500  A.  The  small  sensi¬ 
tive  area  of  the  silicon  cells  nakes  them  unattractive.  Selenium  photo¬ 
cells  are  available  in  a  range  of  sizes  with  quartz  cover  plates.  These 
cells  suffer  from  the  same  shortcomings  of  other  solid-state  photodetec¬ 
tors:  relatively  long  time  constants  and  temperature  coefficients. 
Because  of  the  extreme  cost  and  volume  advantages  of  the  selenium  photo¬ 
cell  over  vacuum  phototubes  and  photomultipliers  it  was  believed  that 
these  should  be  used  if  possible. 
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When  selenium  photocells1  are  loaded  by  a  low- input-impedance 
amplifier,  linear  response  is  obtained  with  low-frequency  signals. 
Response  of  the  cells  to  a  varying  signal  can  be  described  by  two 
characteristic  time  constants,  one  determined  by  the  cell  capacitance 
and  the  load  resistance,  and  the  other  a  fixed  internal  characteristic 
of  the  cell.  The  load  Impedance  can  be  made  low  enough  so  that  the 
first  time  constant  is  small  compared  to  the  sampling  rate.  A  suitable 
compensation  is  applied  to  the  data  to  remove  the  effects  of  the 
internal  time  constant. 

Photometer  Amplifier.  Photocell  signal  currents  are  on  the  order 
of  2  x  10“ 9  amp  and,  as  stated  above,  a  basic  requirement  for  linearity 
is  a  low  load  impedance.  An  operational  amplifier  used  in  a  current - 
to-voltage  converter  mode  is  capable  of  supplying  the  gain  with  low 
input  impedance.  Because  of  the  temperature  range  encountered  during 
the  flight  the  drift  requirements  are  severe.  A  limit  of  A  x  10~12 
amp/*C  is  required  to  limit  output  change  to  0.1  v  over  this  temperature 
range.  This  performance  is  available  from  a  negative  feedback  transis¬ 
tor  amplifier  with  a  capacitive  coupling  to  the  input.  Although  DC 
response  cannot  be  obtained,  with  a  band  pass  from  0.01  to  20  Hz  a 
suitable  compensation  can  be  applied  to  the  data  to  recover  the  pulse 
levels  in  a  quasi-steady-state  situation.  The  method  of  compensation 
is  described  in  the  Appendix. 


ASSOCIATED  PAYLOAD  COMPONENTS 

Telemetry.  The  photometer  output  signal  is  fed  to  a  standard 
subcarrier  oscillator  (voltage-controlled-oscillator)  with  a  10.5  kHz 
center  frequency.  That  output  is  suitably  amplified  and  modulates  the 
transmitter,  which  is  of  the  type  used  in  meteorological  rocket  sounding 
payloads.  The  frequency  is  1,680  mHz.  A  spiral  helix  slot  antenna  is 
used  to  give  a  reasonably  isotropic  radiation  field. 

Package  Temperature.  To  monitor  the  photometer  sensing  head 
temperature  during  flight  a  bead  thermistor  is  installed  behind  the 
photocell.  It  is  part  of  a  shunted  voltage  divider  that  has  an  output 
voltage  in  the  range  of  0  to  5  v.  This  teupe  ature  is  sampled  for 
1.5  sec  every  5  min  by  an  electronic  commutator  in  the  photometer  signal 
circuit. 

Marker  Pulse  Generator.  To  allow  for  automatic  digital  recording 
readout  of  the  telemetered  data,  marker  pulses  are  generated  between 
each  pair  of  light  measurements.  The  system  consists  of  a  micro- 


1  The  selenium  photocells  used  for  this  application  were  manufactured 
by  Dr.  Bruno  Lange  in  Berlin,  Germany. 
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miniature  light  bulb  mounted  in  the  photometer  case  above  the  filter 
wheel.  Small  holes  in  the  filter  wheel  allow  brief  pulses  of  light  to 
pass  to  a  photodiode.  The  photodiode  is  used  as  a  gate  to  allow  bursts 
of  signal  from  a  7-kHz  oscillator  to  pass  to  an  amplifier.  The  amplifier 
output  is  then  mixed  with  the  subcarrier  oscillator  signal. 

Power  Supply.  For  primary  power  four  3-amp-hr-silver-zinc  cells  are 
connected  in  series.  The  electronic  components  require  +28,  -28,  and 
+125  vDC.  These  are  supplied  by  a  DC-to-DC  converter.  The  transmitter 
filament,  filter  wheel  motor,  and  miniature  light  source  are  supplied 
directly  from  the  primary  power  source. 


PHYSICAL  CONSTRUCTION 

Packaging  and  Thermal  Control.  All  components  are  stacked  such  that 
the  maximum  diameter  conforms  to  an  extended  Areas  nose  cone.  For  thermal 
stability,  shock  and  vibration  control,  and  component  protection,  the 
active  assembly,  exclusive  of  photometer  sensing  head,  is  potted  in  a 
medium  density  polyurethane  foam.  To  ensure  dimensional  stability  of  the 
potting  during  upflight  this  foam  is  cured  at  70*C  for  12  to  24  hr.  All 
nonconmerclal  electronic  circuits  are  potted  in  Eccogel  1265  epoxy  resin 
for  thermal  and  mechanical  stability.  Figure  8  is  a  drawing  of  the 
completed  assembly. 

Deployment  System.  To  minimize  the  problems  of  obscuration  of  the 
sun  sad  backscattered  light  from  the  parachute,  a  device  has  been  con¬ 
structed  to  lower  the  instrument  approximately  15  ra  below  the  base  of 
the  parachute.  This  device  consists  of  a  coil  of  3/16-in.  nylon  rope 
installed  in  a  deployment  can  that  mounts  on  the  photometer  head.  At 
ejection  this  can  stays  with  the  parachute  and  the  instrument  falls  away 
until  the  rope  is  completely  extended. 

Ground  Station  and  Digital  Data  Recording.  As  noted  previously  due 
to  the  high  initial  descent  rate  at  altitudes  at  which  data  are  particu¬ 
larly  desirable  a  high  sampling  rate  is  required.  In  the  interest  of 
simplicity  the  same  sampling  rate  is  maintained  throughout  the  flight. 

This  mass  of  data  points  is,  of  course,  useful  but  manual  data  readout 
is  prohibitively  time-consuming.  For  this  reason  an  automatic  digital 
recording  system  has  been  constructed.  The  demodulated  telemetry  signal 
is  taken  from  the  GMD  meteorological  receiver  after  the  detector.  A 
line  driver  has  been  constructed  for  insertion  in  the  socket  of  the 
V-1001  modulation  amplifier  and  cathode  follower.  The  signal  consists 
of  the  frequency-modulated  10.5-kHz  subcarrier  channel  mixed  with  7-kHz 
marker  pulses.  These  frequencies  are  separated  by  band  pass  filters  and 
are  fed  to  a  logic  assembly  (Fig.  9).  The  marker  pulses  generate  gating 
pulses  that  control  two  parallel  counters.  One  counter  accumulates  the 
number  of  crossovers  of  a  i2.5-kHz-crystal-controlled  oscillator  in  the 
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ground  station,  used  as  a  measure  of  elapsed  time  between  marker  pulses. 
The  second  counter  accumulates  the  number  of  subcarrier  crossovers  in 
the  same  interval.  During  the  marker  pulse  the  previous  counts  are 
transferred  to  storage  registers,  the  counters  are  reset  to  zero,  and  a 
new  count  begins.  As  the  new  count  is  accumulating,  the  previous  time 
and  data  words  are  serialized  and  punched  out  on  paper  tape  together 
with  a  parity  channel  and  flags  to  indicate  time  or  data  words.  In 
order  to  provide  a  continuous  time  count  so  that  elapsed  time  from 
launch  can  be  maintained  in  relation  to  real  time  an  automatic  dump  mode 
is  available.  This  is  required  due  to  dropouts  in  the  telemetry  signal 
that  oc  cur  on  the  upflight  and  after  ejection  before  the  parachute- 
payload  system  stabilizes.  In  the  automatic  dump  mode  if  a  marker  pulse 
' 3  not  received  before  the  time  counter  fills  up  to  a  present  level  a 
dump  to  the  tape  punch  is  made  and  an  unique  flag  is  punched. 

Auxiliary  systems  consist  of  a  discriminator  and  strip-chart  recorder 
for  in-flight  monitoring  of  performance  and  an  analog  tape  recorder  on 
which  is  recorded  the  raw  demodulated  telemetry  data,  the  12.5-kHz  timing 
frequency,  a  standard  real-time  code,  and  voice  annotation.  The  discrim¬ 
inator  serves  a  dual  purpose  in  that  the  channel  band-pass  filter  serves 
as  the  data  channel  filter  for  the  digital  system. 

This  type  of  digital  data  recording  system  has  the  advantage  of  high 
frequency  noise  rejection  because  of  integration  over  a  number  of  cycles 
that  average  to  zero.  Disadvantages  are  sensitivity  to  impulse  noise 
and  telemetry  dropouts. 


DATA  ANALYSIS 


The  procedure  for  reduction  of  ROCOZ  flight  data  to  ozone  distribution 
is  lengthy  and  is  the  primary  shortcoming  of  the  optical  technique. 
However,  because  the  method  yields  absolute  ozone  concentrations,  albeit 
with  smoothing  over  any  detailed  structure,  the  task  must  be  considered 
worthwhile. 

A  flow  chart  of  the  complete  analysis  process,  including  calibration 
is  shown  in  Fig.  10.  Although  it  is  possible  to  handle  all  steps  in  a 
computer,  at  this  time  only  a  certain  number  have  been  programmed.  These 
are  the  steps  involving  the  most  data  handling. 

The  flight-data  reduction  process  starts  with  either  conversion  of 
the  punched  paper  tape  generated  in  flight  to  a  computer-compatible 
digital  tape  or  by  digitization  of  the  analog  tape.  This  digital 
magnetic  tape  contains  only  the  crossover  counts  of  the  data  channel 
and  a  reference  constant  frequency  oscillator  between  marker  pulses. 

To  convert  these  counts  to  relative  irradiance  and  elapsed  time  a 
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second  program  is  required.  The  inputs  to  this  program  are  the  data 
tape  and  the  compensation  matrix  derived  from  the  preflight  photometer 
calibration.  Because  this  matrix  is  a  function  of  temperature  the  one 
corresponding  to  the  average  package  temperature  in  flight  is  selected. 
The  outputs  of  this  program  are  the  crossover  counts  per  unit  time  for 
each  filter  wheel  position,  the  compensated  intensity  values  for  the 
four  filters,  and  the  ratios  of  the  absorbed  wavelength  intensities 
(2700,  2950,  and  3050  X)  to  the  reference  wavelength  intensity  (3300  X). 
These  values  are  stored  and  listed  for  each  sample  together  with  the 
elapsed  time  at  the  beginning  of  each  frame  from  the  start  of  the 
recording. 


FIG.  10.  R0C0Z  Data  Processing. 
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The  next  step  In  the  reduction  procedure  is  to  filter  the  data  to 
eliminate  the  noise  contributions  due  to  system  noise  and  pendulation 
on  the  parachute.  To  assess  the  filtering  parameters  as  well  as  to 
evaluate  the  over-all  operation  and  compensation,  the  raw  data  are 
plotted  vs.  time  by  a  computer  program.  Two  filtering  routines  have 
been  used  with  ROCOZ  data.  The  first  was  a  low-pass  digital  filter  in 
which  the  cut-off  frequency  and  sharpness  are  adjustable  parameters. 

Because  the  natural  frequency  of  the  intensity  vs.  height  (or  time) 
curve  is  very  low  and  the  mechanics  of  the  system  introduces  transients 
in  the  data,  this  type  of  filter  was  not  very  successful. 

The  filter  presently  being  used  is  a  least-squares  sliding  polynomial 
curve-fitting  routine.  It  has  the  advantages  of  rapid  response  to 
transients  (which  can  then  be  appropriately  treated)  and  a  true  response 
to  clean  data.  The  output  of  this  filter  is  presently  printed  and  plotted. 
It  then  becomes  necessary  to  decide  whether  the  filtering  parameters  have 
been  selected  correctly  and  whether  additional  smoothing  is  possible  or 
desirable.  In  general  with  the  present  system,  because  of  transients 
introduced  by  the  dynamics  while  on  the  parachute,  it  is  then  necessary 
to  manually  smooth  the  filtered  data  to  remove  the  ficticious  transients. 
The  result  is  intensity  vs.  time  curves  for  each  sampled  wavelength. 

It  is  normally  most  convenient  to  work  with  the  intensity  data  as  a 
function  of  height  rather  than  time.  The  height  vs.  time  data  are 
obtained  independently  from  a  tracking  radar  and  are  typically  furnished 
to  the  user  as  a  computer  printout.  Irregularities  are  removed  by 
smoothing  a  manual  plot  of  the  data.  Times  for  crossing  selected  heights 
(typically  each  kilometer)  are  then  read  off  and  used  to  construct  inten¬ 
sity  vs •  height  profiles. 

Additional  input  data  required  to  calculate  the  ozone  distribution 
are  the  effective  absorption  coefficients  for  each  filter,  sec  z  vs. 
height,  and  air  temperature  vs.  height. 

The  effective  absorption  coefficients,  a j ,  are  calculated  using  Cq.  3 
that  has  been  programmed  for  integration  with  a  wavelength  increment  of 
1  k.  While  this  wavelength  resolution  is  greater  than  is  required,  the 
determination  of  input  parameters  with  fine  structure  (solar  spectrum, 
ozone  absorption  coefficient)  is  made  easier  in  that  the  curve  reader  is 
not  required  to  determine  averages  over  the  wavelength  increment  with 
great  accuracy.  It  is  assumed  that  reading  errors  are  random. 

The  values  of  the  input  parameters  were  obtained  as  follows:  The 
extraterrestrial  solar  spectrum,  IQ(X),  was  read  from  a  tracing  of  data 
obtained  by  the  Haval  Research  Laboratory2  on  14  June  1949.  This  trace 


2  Furnished  to  us  by  Dr.  R.  Tousey. 
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was  read  at  1  A  intervals  in  the  range  2650  to  3399  A.  Because  the 
original  trace  was  presented  in  terms  of  log  I  vs.  X  a  short  program 
was  writcen  to  convert  to  1(A)  and  change  to  the  card  format  used  in 
the  0^  program. 


The  ozone  absorption  coefficient  values  a (A)  were  read  from  an 
expanded  scale  plot  of  Inn  and  Tanaka's  (Ref.  1)  and  Vigroux's  (Ref.  6) 
data  and  punched  on  cards  in  the  necessary  formats.  The  18*C  values 
were  used.  The  Rayleigh  scattering  (extinction)  data  could  not  be 
obtained  in  a  tabulation  suitable  for  this  application  so  a  program  was 
written  for  calculation  from  the  equation; 


_  8ir3  [(GF  +  l)2  -  l]2  D  l<r6 
6<X)  1 - P - F  6.4894  ” 


in  which 


G  -  10-8/6432.8  +  +  -“5^0 

146-121  41  -121 

A2  A2 


,  and 


P[1  +  (0.7870  -  0.0118  T)10~6  P] 
960.857(1  +  0.C03661  T) 


where 


P  »  atmospheric  pressure  in  millibars 
T  -  atmospheric  temperature  in  *C 
D  •  1.061,  the  depolarization  factor  for  air 
N  -  relative  number  density 
M  «■  total  air  mass  in  molecules/cm2 


The  filter  transmission  curves  are  obtained  by  measurement  on  a  Cary 
Model  14  spectrophotometer.  The  transmission  values  at  each  wavelength 
were  originally  read  from  the  tracings  but  presently  are  directly 
obtained  in  digital  form  by  use  of  an  analog  to  digital  converter.  A 
computer  program  has  been  written  to  compensate  for  variations  in  the 
100%  transmission  level  and  shifts  in  the  zero  level  and  to  convert  the 
punched  paper  tape  output  of  this  device  to  cards  in  the  appropriate 
format.  A  similar  technique  is  used  for  input  of  the  sphere  transmission 
data. 


The  photocell  spectral  response  curve  was  measured  in  the  Cary 
spectrophotometer  by  comparison  with  the  fluorescent  light  of  a  sodium 
salicylate  coated  filter  disc.  The  data  being  used  match  with  that, 
furnished  by  Paetzold3  in  the  crossover  region  at  the  long  wavelength 


3  Private  communication. 
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end  and  have  a  form  similar  to  published  data.  As  a  check  on  the  error 
possibly  involved  by  use  of  the  wrong  curve,  the  u^’s  were  calculated 
using  a  flat  spectral  response  over  the  entire  wavelength  region.  The 
results  differed  by  less  than  12  from  those  obtained  with  the  assumed 
response  curve. 

The  output  of  the  cu  program  is  a  listing  of  ctj  and  3aj/3x  for  a 
sequence  of  increments  in  x  and  M  over  the  ozone  and  air  mass  ranges 
suitable  for  e?~h  filter.  These  values  are  plotted  for  ease  of  manual 
ozone  computation. 

The  values  of  sec  z  are  calculated  from  position  and  time  data  with 
reference  to  the  Air  Almanac  for  solar  GHA  and  declination  using  the 
equation 

— - —  •  cos  z  «  cos  LHA  cos  6  cos  6  +  sin  6  sin  6 
sec  z  T  T 

where 

♦  -  latitude 
6  «  solar  declination 
LHA  •  GHA-longitude 

For  use  here  the  values  are  calculated  at  convenient  times,  plotted,  and 
the  required  values  read  from  this  curve.  (More  recently  a  computer 
program  has  been  written  for  this  calculation.) 

From  these  parameters  and  the  intensity  distributions  the  ozone 
distribution  is  calculated  using  Eq.  5.  The  convergence  is  rapid,  such 
that  only  two  iterations  are  required.  It  is  possible  to  make  precise 
corrections  for  the  ozone  absorption  coefficient  temperature  dependence 
in  two  cases:  (1)  when  the  concentration  decreases  rapidly  with  height 
as  above  the  main  peak,  end  (2)  for  Isothermal  regions.  Although  the 
techniques  for  introducing  this  temperature  correction  have  not  been 
formalized,  the  present  procedure  is  to  evaluate  the  distribution  assuming 
18*C  ozone,  then  apply  the  corrections  derived  from  Vigroux's  data.  The 
air  temperatures  used  in  this  calculation  are  obtained  from  independent 
radiosonde  or  rocketsonde  flights. 

The  altitude  increment  size  used  in  the  ozone  calculation  determines 
the  resolution  of  the  distribution.  The  quality  of  the  intensity  data 
is  used  as  a  guide  in  selection  of  the  minimum  increment.  Generally, 

2  km  is  used  unless  data  of  exceptional  quality  are  obtained. 
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FLIGHT  TESTING 


The  rocket-test-flight  program  of  the  complete  system  was  initiated 
in  December  1963  after  several  balloon  flights  of  prototype  instruments 
in  which  optical  filter  parameters  and  the  mechanical  design  were  evalu¬ 
ated.  An  earlier  rocket  flight,  in  June  1963,  was  made  to  test  the 
deployment  system.  With  the  exception  of  a  single  flight  at  White  Sands 
Missile  Range,  the  tests  were  conducted  at  the  Pacific  Missile  Range 
(PMR) ,  Point  Mugu,  California.  All  flights  used  the  system  described 
in  this  report  except  for  minor  changes  in  potting  procedures. 

The  ground  station  for  the  earlier  flights  consisted  basically  of  an 
analog  tape  recorder,  crystal-controlled  oscillator,  discriminator, 
strip-chart  recorder,  and  a  switching  panel.  This  system  was  used  for 
primary  recording  until  the  digital  recorder  became  available  later  on 
in  the  test  series.  At  that  time  its  use  was  continued  for  back-up 
purposes.  Tape  recorder  channel  functions  were  voice  annotation, 
iKIG  B  timing  signal,  telemetry  signal,  and  12.5  kHz  signals.  The 
strip-chart  record  was  used  to  .unitor  in-flight  performance. 

During  the  entire  flight  series  the  analog  tape  records  were 
digitized  using  the  data  automation  system  at  this  Center.  For  conver¬ 
sion  of  the  real-time  punched  paper  tapes  to  computer-compatible  magnetic 
tapes,  a  program  was  written  by  PMR  Weather  Center  personnel  for  use  on 
their  meteorological  computer.  These  tapes  were  then  returned  to  China 
Lake  for  further  reduction. 

A  listing  of  test  flights  and  various  performance  characteristics 
is  given  in  Table  1.  This  program  suffered  greatly  from  rocket  problems 
with  three  motor  malfunctions  out  of  the  11  flights.  After  higher  quality 
rockets  were  obtained  problems  in  the  transmitter  and  payload  potting 
procedures  resulted  in  further  failures.  These  problems  were  resolved 
in  the  final  flight  in  October  1964.  An  analysis  of  this  flight  is 
given  in  the  next  section  of  this  report.  Although  this  flight  was  made 
with  less  than  optimum  solar  zenith  distance  conditions,  it  was  possible 
to  examine  the  system  performance  and  derive  an  ozone  distribution. 

Shortly  after  the  test  flights  were  completed,  a  series  of  flights 
was  made  from  the  U.  S.  S.  Croatan  on  the  NASA  Mobile  Launch  Expedition 
41.  From  these  nine  flights  it  is  possible  to  make  a  realistic  estimate 
of  the  package  reliability  and  set  certain  conditions  for  its  use.  Of 
these  flights,  one  rocket  motor  and  one  ground  station  (GMD)  failure 
occurred.  Of  the  payloads,  one  transmitter  failed  due  to  over-heating 
during  an  extended  prelaunch  hold  and  an  integrating  sphere  was  broken 
in  another  flight.  The  net  result  was  five  successful  flights  out  of 
nine  launches. 


TABLE  1.  ROCOZ  Test  Flight  Suanary 
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ANALYSIS  OF  DATA  FROM  FLIGHT  45 


The  final  flight  In  the  test  series,  NWC  Flight  45,  was  launched  at 
PMR  on  8  October  1964  at  17b  40«  58s  GMT  with  a  solar  elevation  angle 
of  41  degrees.  The  flight  operations  were  conducted  by  the  Geophysics 
Division  of  PMR.  Peak  altitude  achieved  was  163,000  ft  (49.8  km)  which, 
for  the  payload  weight,  is  somewhat  lower  than  nominal  180,000  ft  (55  km). 
Figure  11  shows  the  altitude  vs.  time  for  this  flight.  Optical  signals 
from  the  photometer  were  received  within  0.2  sec  after  ejection  that 
took  place  at  1743  hr.  The  parachute  opened  1.4  sec  after  ejection  and 
a  clean  photometer  signal  was  received.  This  signal  during  the  ejection 
and  deployment  sequence  is  shown  in  Fig.  12.  The  data  frames,  consisting 
of  pulses  from  the  four  filters  and  a  zero  level,  showed  variability  in 
relative  pulse  height  due  to  the  rapid  rotation  of  the  payload  soon  after 
ejection.  Within  10  sec  after  ejection  a  longer  period  modulation  of 
frame  amplitudes  also  appeared.  This  modulation  was  regular  with  an 
initial  period  of  8  sec  and  amplitudes  up  to  80%.  The  modulation  contin¬ 
ued  with  high  amplitude  until  1811  hr.  Because  of  the  low  solar  elevation 
angle,  pendulation  angles  exceeding  11  decrees  would  result  in  modulation 
of  the  type  observed.  An  increase  in  amplitude  of  the  3050-  and  3300-A 
pulses  occurred  at  1750  hr  such  that  clipping  of  the  highest  pulses  took 
place  until  1802  hr.  Telemetry  was  generally  clean  throughout  the  flight 
although  large  fluctuations  in  the  RF  signal  strength  were  observed  ini¬ 
tially,  decreasing  as  the  package  descended.  The  data  were  recorded  on 
analog  magnetic  tape  and,  using  the  logic  system,  on  punched  paper  tape 
in  digital  form.  Because  of  the  signal  strength  variations  some  diffi¬ 
culty  was  noticed  in  maintaining  the  proper  level  for  triggering  of  the 
logic  system  by  the  marker  pulses. 

The  package  temperatures  were  read  out  from  the  in-flight  analog 
strip-chart  record  and  are  shown  in  Fig.  13.  A  27 *C  total  change  was 
observed,  but  during  the  portion  of  the  flight  in  the  primary  ozone 
region  above  20  km  the  changes  are  limited  to  11*C.  For  this  reason 
the  compensation  matrix  for  room  temperature  was  used  for  the  data 
readout . 

Two  methods  of  digital  data  readout  were  used  to  provide  a  comparison 
of  real-time  analog  and  digital  recording.  The  analog  magnetic  tapes 
were  processed  by  the  data  automation  system  at  this  Center.  The  punched 
paper  tapes  were  read  and  processed  in  part  by  the  PMR  Geophysics  Division 
computer  personnel.  No  useable  data  could  be  obtained  from  either  source 
prior  to  1800  hr  due  to  excessive  modulation  of  the  optical  signal  and 
loss  of  marker  pulses  during  weak  portions  of  the  telemetry  signal.  The 
data  from  the  analog  magnetic  tape  were  significantly  lower  in  quality 
than  that  from  the  punched  paper  tape.  From  an  examination  of  the  analog 
strip-chart  records  it  appears  that  useable  data  can  be  obtained  from  the 
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real-time  digital  records  when  the  modulation  of  the  light  on  the 
photometer  is  less  than  30%  in  amplitude  (i.e.,  the  steady-state  approx¬ 
imation  holds  under  those  conditions).  For  greater  modulation  levels 
the  compensation  method  diverges. 

To  obtain  intensity  vs.  height  data  at  higher  altitudes  the 
amplitudes  of  selected  pulses  on  the  strip-chart  record  were  measured 
manually.  Smooth  curves  have  beta  fit  to  these  points.  The  results 
are  given  in  Fig.  14  where  the  light  intensities  at  the  three  shortest 
wavelengths  relative  to  the  reference  wavelength  are  shown.  The 
clipping  problem  resulted  in  loss  of  data  in  the  20  to  30  km  region. 

The  characteristics  of  the  optical  filters  used  for  this  flight  are 
given  in  Table  2.  The  listed  parameters  are  center  wavelength,  Ac, 

width  at  the  50%  relative  transmission  poirsts  AA^.  peak  transmission 
Tjnax*  effective  band  pass  AA j ,  observed  and  calculated  extraterrestrial 
signal  from  the  photometer  relative  to  the  3300  X  signal  S^'S  .  The 
discrepancy  between  observed  and  calculated  ratios  s*s  not  been  explained. 


TABLE  2.  Optical  Filter  Parameters 


j  Filter  No. 

V 

A 

AA^, 

X 

T 

max 

AA 

A 

S,/S„  obsd. 

J  ° 

S./S  calcd. 

j  ° 

0 

3287 

47 

0,0041 

0.193 

•  ♦  ♦  • 

•  •  •  • 

1 

3100 

48 

0.0267 

1.28 

1.18 

3.93 

2 

2973 

54 

0.0301 

1.62 

0.97  (est.) 

3.41 

3 

2741 

] 

50 

1 

0.0738 

3.70 

•  •  •  • 

2.05 

The  effective  ozone  absorption  coefficients  vs.  total  ozone  amount 
for  this  filter  set  are  shown  in  Fig.  15.  The  slight  influence  due  to 
Rayleigh  scattering  is  shown  by  the  curves  for  different  values  of  m, 
the  relative  air  mass  above  the  instrument.  Here,  only  the  influence 
of  single  scattering  is  taken  into  account.  The  absorption  coefficient# 
are  computed  for  ozone  at  18*C.  Especially  in  the  lower  stratosphere 
the  ozone  temperature,  assumed  in  equilibrium  with  the  air,  is  signifi¬ 
cantly  lower  than  18*C.  The  temperature  dependence  of  the  ozone  absorp¬ 
tion  coefficient  is  on  the  average  a  slowly  changing  function  of  the 
wavelength  and  a  multiplicative  temperature  function  can  be  used  for 
each  filter. 

The  ozone  concentrations  were  calculated  using  Eq.  5  and  the 
associated  slope  correction  (E-].  6).  The  value  of  sec  z  ranged  from 
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FIG.  15. 


Curves. 


1.527  at  launch  to  1.456  at  22  km.  Considering  the  data  quality  an 
altitude  increment  of  2  km  was  selected.  Evaluation  was  not  attempted 
below  30  km  because  of  the  previously  mentioned  clipping  problem. 

Results  are  shown  in  Fig.  16.  Also  plotted  are  the  results  of  an  earlier 
measurement  at  White  Sands  Missile  Range  on  14  June  1949  as  corrected  by 
Craig  (Ref.  8).  A  detailed  error  analysis  has  not  been  conducted  at  this 
time;  however,  an  estimate  can  be  made  from  the  concentrations  derived 
from  two  filters  at  the  crossover  region.  In  this  case,  at  43  km  a  22% 
difference  exists.  It  is  expected  that  better  agreement  will  be  obtained 
under  more  optimum  flight  conditions. 

The  results  of  this  flight  demonstrate  that  (1)  the  basic  system  can 
operate  successfully,  (2)  the  payload-parachute  system  must  be  stabilized 
to  acquire  digital  data  and  higher  quality  analog  data  over  the  entire 
altitude  range,  and  (3)  the  real-time  digital  recording  system  produces 
higher  quality  data  than  a  digitized  analog  magnetic  tape  system. 
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Appendix 

LINEAR  OVERSHOOT  MODEL  OF  AMPLITUDE  COMPENSATION 


In  the  ROCOZ  system  the  irradiance  in  four  narrow  spectral  bands 
and  a  zero  light  level  are  sampled  sequentially.  Each  group  of  five 
samples  is  denoted  a  frame  and  for  steady-state  illumination  conditions 
the  frames  are  identical  and  occur  repetitiously.  The  time  variation 
of  the  irradiance  in  each  pulse  at  the  photodetector  is  determined  by 
the  system  geometry  and  consists  of  a  symmetrical  pulse  with  a  rounded 
top.  In  a  direct-coupled  system  with  a  fast-response  detector  the 
individual  pulse  amplitudes  or  the  pulse  areas  serve  as  a  measure  of 
the  Irradiance  when  referred  to  the  zero  light  level  value. 

In  this  application,  due  to  cost  and  physical  space  considerations, 
the  use  of  a  solid-state  photodetector  is  desirable.  This  device  has 
internal  time  constants  of  the  same  order  of  magnitude  as  the  sampling 
period  that  is  determined  by  flight  parameters.  A  further  complication 
arises  from  the  low  irradiance  levels  available.  Due  to  the  small 
signal  currents,  suitable  direct-coupled  amplifiers  with  satisfactory 
drift  specifications  have  not  been  available  until  recently. 

It  has,  therefore,  been  necessary  to  derive  a  method  for  compensation 
of  the  observed  signals  to  obtain  the  input  irradiances.  Although  the 
method  is  illustrated  with  analog  wave  forms,  it  is  applied  to  a  discrete 
set  of  measurements  that  consist  of  the  areas  under  the  pulses.  Each 
pulse  is  characterized  by  its  position,  i,  in  the  frame.  Residual 
effects  persisting  after  the  input  pulse  ha»  gone  to  zero  are 
characterized  by  the  integer  j  equal  to  nAi. 

An  input  pulse  of  light,  1^  is  transformed  by  the  photocell  and 
amplifier  to  the  output  pulse 

k,i,j  k,i,j  k,i 

where  Aj^i  is  the  average  signal  level  over  one  or  more  frames  preceding 
the  sample  Ij^i  (Fig.  17). 
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FIG.  17.  Amplifier  Overshoot  Characteristics. 


Each  output  pulse  is  modified  by  the  preceding  pulses  such  that 

Sk,i,0  "  BIk,i  "  Sk,i-l,l  ”  Sk,i~2,2  ”  Sk,i-3,3 
n 

“  BIk,i  '  S  Sk,i-j,j 

j-1 

where  if  j>i,  set  i-j-*i-j+5  and  k*k-i,  and  B  is  defined  as  a  sensitivity 
constant.  In  a  steady-state  situation,  it  is  obvious  that  Sk(ij 
for  j>4  for  a  five  sample  per  frame  device.  Thus  in  this  case  there  is 
the  following  set  of  equations 
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DETERMINATION  OF  THE  OVERSHOOT  COMPENSATION  PARAMETERS 


The  compensation  parameters  (j  t  0)  could,  in  principle,  be 

determined  from  the  transfer  function  of  the  photometer4  but  in  practice 
an  empirical  approach  has  proved  useful. 

Five  input  light  pulses  of  known  amplitude  are  fed  to  the  photometer. 
Consider  the  simple  case  where  1^  *  1^  *  I^  ■  I^  *  0  and  1^  **  I.  The 

subscript  k  has  been  deleted  because  in  a  steady-state  situation  all 
frames  are  identical.  The  value  of  I  is  set  at  various  values  throughout 
the  useful  range  of  operation. 


It  is  assumed  that  the  overshoot  amplitude  depends  only  on  the 
amplitude  of  the  initiating  pulse  and  that  the  overshoot  from  a  zero 
amplitude  pulse  is  zero.  The  five  equations  reduce  to 


S0,0  “  _SA,1 
Sl,0  "  "S4,2 
S2,0  "  ~S4,3 

S3,0  “  *S4,4 

S.  _  -  BI. 

4,0  4 

If  it  is  also  assumed  that  the  overshoot  values  do  not  depend  on  the 
position  of  the  originating  pulse  in  the  frame,  then 

S0,0  “  ”S1  =  '&1 

Sl,0  “  -S2  5  -&2 

S2,0  *  'S3  2  ‘*3 

S3,0  "  "S4  2  '&4 


S4,0  *  81 


4  R.  J.  Stirton,  Test  Department  Technical  Note  3038-55,  January  1964, 
limited  distribution. 
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The  functional  dependence  Aj(I)  is  obtained  experimentally. 

Figure  18  shows  a  typical  response.  It  should  be  noted  a  linear  system 
is  required.  The  output  pulse  amplitude  S^  ^  vs.  I  is  shown  in  Fig.  19. 
This  response  is  also  linear.  * 


-Aj  *Kj  s 


234  3*76*0 

RELATIVE  OUTPUT  SI6NAL 


FIG.  18.  Overshoot  Compensation  Parameters  vs.  Signal  Level. 


INPUT  SIGNAL  LEVEL 


FIG.  19.  Signal  Amplifier  Linearity. 
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The  compensation  functions  can  be  written  as 

yv  •  Vi 

where  the  value  of  K.  is  determined  from  a  calibration,  and  the  full  set 


of  equations 

becomes 

l 

S„  - 

am 

BI  _  - 

K,  I, 

-  K-I- 

-  K-I- 

-  K.  I, 

0,0 

0 

14 

2  3 

3  2 

4  1 

S,  „ 

BI,  - 

K,  I„ 

-  K_I. 

-  I- 

-  K,  I- 

1,0 

1 

10 

2  4 

3  3 

4  2 

S„  „ 

m 

BI.  - 

K,  I, 

-  K-I- 

-  K.  I. 

-  K.  I. 

2,0 

2 

1  1 

2  0 

3  4 

4  3 

S„  „ 

m 

BI_  - 

K,  I„ 

-  K_I, 

-  K-I- 

-  K.  I. 

3,0 

3 

12 

2  1 

3  0 

4  4 

S,  „ 

BI,  - 

K,  I- 

-  K-I- 

-  K-I, 

-  K,  T 

4,0 

4 

13 

2  2 

3  1 

4  0 

In  the  present  application  Iq  ■  0  and  when  rearranged  the  set 
reduces  to  the  nonhomogeneous  set  of  linear  equations: 

Vl  +  K3x2  +  K2x3  +  K1X4  -  "S0,0 

-BI!  +  K4X2  +  K3X3  +  K2X4  “  "Sl,0 

K1X1  -  BI2  +  K4X3  +  K3X4  “  -S2,0 


Vl  +  K1X2  '  BI3  +  K4X4  *  -S3,0 


K3!1  v  K2I2  +  -  BI4  -  -S 


4,0 


where  the  first  equation  can  be  considered  redundant.  For  simplicity, 
it  is  assumed  that  B  ■  1.  In  matrix  form, 

KI  -  -S„ 
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fe 

Is  where 

!*■ 

| 


£ 

■-1  K4  K3  k2- 

V 

Sl,0 

f 

K,  -1  K,  K, 

I„ 

s,  . 

~ 

1  4  3 

2 

2,0 

f  K  - 

,  I  “ 

.  -s0  ■ 

k2  Kx  -I  Ka 

X3 

S3,0 

; 

K  K,  -1 

I. 

s,  „ 

3  2  1  J 

L  4J 

L  4»°J 

To  epply  the  compensation  to  general  steady-state  data,  the 
following  is  required: 


I  -  K 


-1 


i 

?_ 


* 

? 


This  compensation  analysis  is  quite  general  with  the  limitation  that 
the  system  be  linear.  Any  number  of  pulses  per  frame  is  allowable 
although  for  practical  purposes,  inversion  of  matrices  greater  than 
7x7  becomes  prohibitive. 


As  a  simple  example  of  the  use  of  this  method,  consider  a  direct- 
coupled  system  with  fast  response.  In  this  case  “  constant  for  all  1 
and  for  the  ROCOZ  system  with  four  samples  and  a  zero  measurement  per 
frame,  ■  1/4.  A  particularly  simple  result  is  obtained: 


j  *  i 


For  the  ROCOZ  photometer,  it  has  been  impossible  to  define  a  unique 
set  of  Kj/s  because  of  variability  in  the  photocells  and  amplifiers 
and  consequently  Individual  calibrations  are  required. 
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